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Abstract 

Among R 3 -valued triples of random vectors (X, Y, Z) having fixed 
marginal probability laws, what is the best way to jointly draw (X, Y, Z) 
in such a way that the simplex generated by (X, Y, Z) has maximal 
average volume? Motivated by this simple question, we study opti- 
mal transportation problems with several marginals when the objective 
function is the determinant or its absolute value. 

Keywords: Optimal transportation, multi-marginals problems, deter- 
minant, disintegrations. 

1 Introduction 

Given two probability measures \i\ and \i2 on M. d and some objective function 
H : R d x M. d — > R, the classical Monge-Kantorovich optimal transportation 
problem consists in finding a probability measure 7 on x M. d having \i\ 
and \i2 as marginals (i.e. a transportation plan between \x\ and ^2) maxi- 
mizing the total objective J RdxK d H(x,y)d'y(x,y). In his famous article PQ, 
Brenier solved the case H(x,y) = (x,y) and proved (under mild regularity 
assumptions) that there is a unique optimal transportation plan which is 
further characterized by the property of being supported by the graph of 
the gradient of some convex function. Brenier's seminal results have been 
extended to the case of more general costs (see McCann and Gangbo [I]) and 
the subject has received a lot of attention in the last 15 years because of its 
numerous applications in fluid mechanics, probability and statistics, PDE's, 
shape optimization, mathematical economics.... The literature on this very 
active field of research is too vast to give an exhaustive bibliography here, 
we rather refer to the books of Villani [B] and Rachev and Riischendorf [5] 
and the references therein. 

In the present article, we are interested in an optimal transportation 
problem with several marginals. Given d probability measures fix, fid 
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on M rf and an objective function H : ~R dxd — > R, the problem is to find a 
probability measure 7 on M. dxd having fj,i, fj,^ as marginals maximizing 
J K dxd H(xi, Xd)d r y(xi, Xd)- In contrast with the case of two marginals, 
there are few results on the optimal transportation problem when more than 
two marginals are involved with the exception of the article of Gangbo and 
Swi§ch [3] who fully solved the case H(x\, ...,Xd) := Ylij ( x i> x j)- F° r this 
particular problem, Gangbo and Swiech proved existence and uniqueness of 
an optimal transportation plan which is supported by the graph of some 
transport map. In the present paper, we will pay attention here to different 
objective functions, namely: H(x) = det(x) or H(x) = |det(x)|. 

The choice of such functions of the determinant is motivated by the 
following simple question: among random M 3 -valued vectors X,Y, and Z, 
with fixed marginal probability laws, what is the best way to draw jointly 
(X, Y, Z) so that the simplex with vertices (0, X, Y, Z) has maximal average 
volume? Denoting by \ii and ^ the (fixed) probability laws of the 
random vectors (X, Y, Z) , the previous problem amounts to maximize 



1x3 



det(x,y,z)\dr/(x,y,z) 



among joint probability laws 7 having fj,\, Hi and ^3 as marginals. In some 
cases, we will see that solving the problem above amounts to solve the 
simpler problem where | det | is replaced by det and we will study this case 
in dimension d > 2. If d = 2, we have det(x,y) = (x,Ry) (with R the 
rotation of angle — 7r/2), so that, up to the change of variables y' = Ry, the 
optimal transportation problem with the determinant is a special case of the 
problem solved by Brenier [I]. 

Let us define some notations. In the sequel, given X a locally compact 
separable metric space, we denote by M.(X) (respectively Ai\{X)) the set 
of Radon measures (respectively of Radon probability measures) on X. If 
X and Y are locally compact separable metric spaces, fj, € A4\(X), and 
/ : X — > Y is a Borel map we shall denote by f%\x the push forward of fi 
through / i.e. the element of M\(Y) defined by f$n{B) = n(f~ l (B)) for 
every Borel subset B of Y. If 7 G M(R dxd ) and i G {l,...,d}, 7^7 G M{R d ) 
is called the i-th marginal of 7 (where 7Tj is the i-th canonical projection). 
One can also define 7Tj$7 by: 

f{xi)d(iri$j)(xi) = / f(x i )dj(x 1 ,...,x d ), 

JR d >< d 

for every bounded and continuous function / on M d . Given d probability 
measures on M. d , /zi,...,/^ , we denote by H(fJ>i, Hd) the set of proba- 
bility measures on M. dxd having /ii,...,/^ as marginals. In other words, 
7 G M\(R dxd ) belongs to U(m, /x d ) if and only if 



/ 



f(x i )d>y(x 1 ,...,Xd)= I f(xi)d[ii(xi),Vi = l,...,d, 

Xd 
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for every bounded and continuous function / on R d . Given H £ C°(M. dxd , R), 
the Monge-Kantorovich optimal transportation problem with marginals fj,i, fj, d 
and objective function H then reads as: 

sup / H(x 1 ,...,x d )dj(xi,...,x d ). 

7 en( w ,...,/i d ) JR d >< d 

We shall focus here on the two special cases: 



(MIC) sup / det(x 1 ,...,x d )d>y(xi,...,x d ) 
7en( w ,...,M d ) J^ dxd 

and 

(7W/C) a sup / \det(xi,...,x d )\dj(xi,...,x d ). 

7en(/Lti,...,/i d ) JM dxd 

The paper is organized as follows. In section 2, we solve a particular 
example which actually gives insight on the general case. Section 3 is devoted 
to existence, duality and characterization of minimizers for (MIC). In section 
4, we construct minimizers in the case of radially symmetric marginals. In 
section 5, we give conditions ensuring that (MIC) and (MK) a are in fact 
equivalent. Section 6 contains various remarks regarding uniqueness issues. 



2 An elementary example 

In this section, we study a simple but illustrative example, which actually 
contains most of the ideas necessary for the understanding of the general 
case. Let us consider the problem (MIC) in dimension 3, with [i\ = \i2 = 
//3 = £g, where £g stands for the uniform probability measure on the unit 
ball B of M 3 , 

sup (/ det(x,y,z)dj(x,y,z) ) . (1) 
7 en(£|,£|,£|) \Jb3 J 

Then, the following result holds. 

Theorem 1 The problem ([!]) admits a solution 7 given by, for every con- 
tinuous f : B 3 — > R, 

/ fd ^ = W\l\l f( x >\ x \y> x ^y)^^) dx > ( 2 ) 

J I I J B V J J 

where 

§(x) = {yeS 2 s.t. (x,y} = 0}. 
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Remark 1. Let us remark that, the support of the optimal measure 7 is 
the set of all the triples (x,y,z) G B 3 , such that \x\ = \y\ = \z\ and (x,y,z) 
is a direct orthogonal basis of M 3 . As we shall see later on, this last property 
comes from the fact that the measures are radially symmetric. In addition, 
the definition of 7 through its successive disintegrations has the following 
probabilistic interpretation in terms of conditional laws. Consider 7 as the 
law of a vector (X, Y, Z) of random vectors in M 3 , each of them following an 
uniform law on the unit ball. Then, (|2j) means that the conditional proba- 
bility of Y given X is uniform on S(X) and that the conditional probability 
of Z given (X, Y) is the Dirac mass at ~n^r ■ 

Remark 2. Let us say a word about the case where the objective function 
is H(x,y,z) = I det(x, y, z)\ (volume maximization). For this H, 7 is still a 
maximizer, but we could have chosen as well 

XAY 

for the third variable, and the solution is not unique. Actually, it is not the 
only source of non uniqueness, and we shall discuss this point in the last 
section. 

Proof of theorem [TJ, 

Assuming that 7 is admissible in (pQ), we only need to prove that it is 
optimal. To do so, consider the variational problem 

^w\{S B ^ x)dx+ L 

with £ the set of triples ((p,ip,x) G C°(B,R) 3 such that 

<p(x) + ip(y) + xO) > det(x, y, z), V(x, y, z) £ B 3 . 
As we shall see in next section, ([3]) is dual to ([I]) in some sense. For all 

(ip, V>, x) e £ , and 7 e n(£ 3 , £ 3 , £| ), 

det(x,y,z)dj(x,y,z) < / (<p(x) + ip(y) + x(z)) d^/(x, y,z), 
B 3 Jb 3 

- Jb 1 * f vW*^ + f ^^ dy + / x(z)dz 

and it immediately follows that 



w . , 1>{v)dy+ / X(z)dz), (3) 

{¥>,ipM££ \Jb Jb Jb 



b jb jb 



sup / det(x,y,z)dj(x,y,z) 
n(£|,£|,£3)7s3 



< inf ■— -r ( / (f(x)dx + / 



+ / X(z)dz ) . (4) 
b Jb 
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Now, consider (^Oj^OjXo) given by 

I X I 

Vx G -B, ^o(x) = ^q(x) = xo(z) = — ■ 

Thanks to the Young inequality, we have 

\/(x,y,z) G B 3 , ipo(x) +ip (y) + Xo(z) > \x\\y\\z\ > det(x,y,z), (5) 

so that (c^o, ipo, Xo) €E £. In addition, according to remark [fl if (x,y,z) G 
supp(7), then all the inequalities in ([5]) become equalities and then, 

V(x, y, z) G supp(7), tpo(x) + ^ (y) + Xo(X) = det(x, y, z). (6) 

It implies that (j3|) is actually an equality and 7 is optimal in (pQ) . To complete 
the proof, it just remains to show that 7 is admissible, which is done in 
proposition [TJ 

Proposition 1 Let 7 defined as in theorem B T/ien, 7 G II(£|, £|, £|) . 

Proof. The fact that 717 (fy = £g is obvious by ([2]). Now, let us suppose 
that we proved 7^7 = £g. Then, we can easily deduce the result for the 
third marginal. Indeed, since for every fixed x G B, the map y ^ is 
one-to-one from §(x) to itself and is simply a rotation with angle vr/2, we 
have that for all continuous / : B — > R, 




/(My) 



S(x) 




S(«) 

We then prove that 7^7 = £g. Let / : B — > M be a continuous function. 
Then 



/ M^tm = ml (l /(My) 

\JS(x) 



f{ry) ^_^± ] d u\ a ) ) dr. 

Applying lemma [2] proved in Appendix A, we get 

g fdfaffl = ^J 1 (f^f{ry)(f^ d !^\d7{\y)\dr 






which ends the proof. 



□ 
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Let us notice that the condition ([6]) completely characterizes the solutions 
of the dual problem ([3|) . In fact (see subsection I3.3P , up to the addition of 
constants that sum to 0, {(po,ipo, xo) is the unique solution of Q. Yet, there 
are infinitely many solutions to ([I]), indeed any 7 G n(£|, C#, £g) having 
its support in all the triples {x,y,z) G B 3 , such that |x| = \y\ = \z\ and 
{x, y, z) is a direct orthogonal basis of R 3 is optimal for ([3]) and we claim 
that there are infinitely many such probability measures (see section . 



3 Duality, existence and characterization 

Given d probability measures on I 11 , , we consider the problem 

(MJC) sup / det(xi, ...,x d )dj{x 1 , ...,x d ). 

In this case, a natural assumption on the marginals fi±, fj, d , is the existence 
of {pi, ...,pd) G (1) +oo) d such that: 

\T- y 1 f \x\ Pi 

y — = 1, and a := } / duAx) < +00. (7) 

Defining for all x = (xi, ....,Xd) G R dxd : 

Ho(x) := J2 



d 1 



,=1 ^ 



and using the fact that | det(x)| < Hq{x), we immediately get 

-a < I det{xi, ...,x d )d'y{xi,...,x d ) < a, V7 G n(/zi,..., (j, d ). 

Similarly, defining for all x = (x\, ....,x d ) G M. dxd : 

H{x) := det(x) + H (x), H{x) := det(x) - H {x), (8) 
we remark that for all 7 G n(/ix, 

/ det(x)d7(x) = / Hdj — a = / Hd^y + a. 

The previous remark implies that when if = det, one may without loss of 
generality replace H = det with H > or H_ < in (A4/C). 

A key point in Monge-Kantorovich theory, is to remark that (in a sense 
that will be made precise later), {MIC) is dual to : 



d . 

( p ) , inf , „ yz / 
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where £ is the set of <i-uples of lower-semi continuous functions ((pi, <p d ) 
from R d to R U {+00} such that 

d 

Y / <fi{x i )>det{x 1 ,...,x d ), V ( Xl ,...,x d ) £R dxd . (9) 
i=i 

Let us note that one obviously has inf(D) > sup(.M/C). 
3.1 The compact case 

In this paragraph, for further use, we consider the case of compactly sup- 
ported marginals and of an arbitrary continuous objective function H. De- 
noting by B r the closed ball in M. d , with center and radius r, we assume 
that fix, are supported in B r for some given r > and that H is an 
arbitrary continuous function on B d . We then consider the optimal trans- 
portation problem: 

sup / H(x 1 , ...,x d )dj(xi, ...,x d ). (10) 

7 6n( A t 1 ,...,/Lt d ) JR d * d 

We define its dual by: 

d „ 

, Yl / VidlH (11) 

(ipi,...,(p d )e£r r~: Jb,. 

where f r is the set of ci-uples of continuous functions (ip\, <^<i) from i? r to 
R such that 

^2<Pi(xi) > H(xi,...,x d ), V (xi,...,x d ) G £ r d . (12) 
i=l 

Proposition 2 Assume that fxi, ■■■■■> Hd are supported in the ball B r and that 
H is continuous on B d then both i!0\) and ( liij) admit solutions and 

max (f77J = min PIP. 



Moreover, Ul\) admits a solution (<pi, <pd) such that for all i = l,...,d 
and all Xi S B r , one has: 



ipi(Xi) = sup 



H(xi,...,x d ) - ^(Pjixj) f (13) 



and, for all x € B r : 



mini? < <Pd(%) < max if, < ipdx) < max if — min if, % = 1, ...,cf — 1. 

B§ Bi Bi Bi 



(14) 
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Proof. 

Step 1 : convex duality 

Equip E := C°(B r ,R) d with the sup norm, and define (the linear contin- 
uous operator) A by A((^i, ...,ip d ))(xi,...,x d ) := Ya=i Viixi) for all (931, ...,ip d ) G 
C°(B r ,R) d and all (x u ...,x d ) G B d . Remark now that (jlip can be rewritten 
as: 

inf F(^) + G(AM) (15) 
with F(ip) = J2i=i J Br and ' for a11 ^ e C°(B*,R): 



Gty) 



if V > -ff 
+00 otherwise . 



The dual problem in the usual sense of convex analysis of (|15p (see [2]) is 
then 

sup -F*(A*7)-G*(- 7 ). (16) 

First, we remark that A* (7) = (7riJJ7, ...,^^7). Elementary computations 
then yield: 



F*(A*( 7 )) 



if 7rj)J7 = m for i = 1, d, 
+00 otherwise . 



and 



G*(- 7 ) 



- | M dxd #^7 if 7 > 
+00 otherwise . 



Hence problem (I16D is exactly (I10p . The Fenchel-Rockafellar duality theorem 
(see [2]) implies then that (|1U|) admits solutions and that 

max (TTUD = inf (HU)- 

Step 2: convexification trick 

It remains to prove that the infimum is attained in (|lip . Let (ipi, <^) € 
£ r and define for all x\ € B r : 



^i(xi) = sup I H(xx,...,x d ) -^2<fj(xj) 

{x 2 ,-.,x d )eB^- 1 [ j =2 



(17) 



by construction ipx — Pi an d (^l? ^d) G £ r . Construct then induc- 

tively -02, tpd-l by setting for i = 2, d — 1 and Xj G B r 



sup 



3>< 
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and finally 



sup 



d-l 



H(xi, ...,x d ) - ^2^j{xj) } . 

3=1 



By construction, tpi < ifi and (-01, ipd) € S r : (tpi, ipd) is therefore 
an improvement of (yi, Vd) in problem (jll|) . On the one hand, since 
V'j < Vj) one nas f° r an * = 1) •••! d and Xi G B r 



sup 



F(a;i,...,a; d ) -^^(^) f • 



On the other hand, the converse inequality holds because (^i, ■ ■■■,ip d ) G £r- 
Step 3: existence of minimizers for All]) 

Using the convexification trick of the previous step, we can find a mini- 
mizing sequence of (|11|) . ip n := (<Pi, ...,(p"j) such that for all n, all i and all 
Xi G B r , one has: 



sup 



H(xi,...,x d )-^2ip](a 



(18) 



Noting that the objective in (fTTI) is unchanged when changing (ipi, cpa) 
into (tpi + ai, ^d + c^d) for constants cc, that sum to 0, we may also assume 
that 

min^™ = 0, Vn G N, Vi = l,...,d-l. 
Together with (fTBj) we deduce that 



and 



minff < <zjj < maxff on i3 r , Vn G 



< max — min H. 



Denoting by to the modulus of continuity of H on Bf, we also deduce from 
(|18p . that for every (x,y) G B%, for every n and i one has: 

|^(x)-^(y)l<^(k-y|)- 

Thus, the sequence <p n is bounded and uniformly equicontinuous hence by 
Ascoli's theorem admits some convergent subsequence. Denoting by <p = 
((/?!, ...,ipd) the limit of this subsequence, it is easy to check that (p belongs 



to £ r , solves (fTTj) and satisfies (|T3|) and (|Mj) . 



□ 
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3.2 The general case 

We now go back to (MIC) (i.e. H = det) for general marginals that only 
satisfy (|7|). By suitable truncation arguments, we have the following result, 
which is proved in Appendix B: 

Theorem 2 Assume that (f?|) is satisfied, then both (MK) and (V) admit 
solutions and 

max (MK) = min (D). 

Moreover, (T>) admits a solution (ip\, (fid) such that for all i = l,...,d 
and all Xi G W 1 , one has: 



(pi{Xi) = sup 

(:r : ,) j7 ^eIR dx ( d - 1 ) 



det(xi,...,x d )-^2(pj(xj)f. (19) 



3.3 Extremality conditions 

This paragraph is devoted to optimality conditions for (MK) that can be 
derived from Theorem [2] (again we are in the case H = det here). For 
(yi, ....;y d -i) G (M^- 1 we denote by Af=i Vi tlie vector yi A ... A y d -\ and 
recall that is characterized by the identity: 



det(yi, ....,y d ,x) = 




For (x\, ....,Xd) € (M d ) d and i G {l,...,d}, f\^Xj is defined in a similar 
way. 

At this point, it is useful to remark that the family of l.s.c. functions 
(ipi, (pa) belongs to £ if and only if for every i E {1, d] one has 

> ^(("l) m A B i). V^i,-,^) G (R d ) d . (20) 

It is also obvious that ((pi, ipd) belongs to £ if and only, it satisfies (j20|) 
for some i G {l,...,d}. Let us finally recall that the convexification trick 
ensures that one can always improve an element of £ in the dual problem 
by replacing it by another element of £ that satisfies f)19[) . Hence solutions 
of (T>) have to agree ®f = ^m almost everywhere with potentials that satisfy 
(fT9j) . Obviously if ((pi, ipd) satisfies (fT9|) . each <pi is a convex potential (as 
a supremum of a family of affine functions). With no loss of generality, we 
may therefore restrict ourselves to the subset £ c consisting of the elements 
(ipi, ....,<p<j) G £ such that each potential ipi is convex on M. d . 

By definition of £ and the duality result of theorem^ we deduce that 7 G 
/id) (respectively ((pi,...<pd) G £) solves (MK.) (respectively solves 
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(£>)) if and only if there exists (c^i, ...ipd) G £ (respectively 7 G n(/xi, ...,/zd)) 
such that 

d 

yjyj(^j) = det(xi, 7-a.e. (21) 

i=i 

For all $ := (991, ifd) G £ , let us define 

d 

Kq, := {( Xl ,...,x d ) G (R d ) d : ^(Pj{xj) = det(xi,...,x d )} 

3=1 

and remark that the (possibly empty) set K§ is closed since it is the minimal 
set of some l.s.c. function. Hence we deduce that (|21 j) is equivalent to 
7 having its support included in K§. If $ G £ c , we have the following 
characterization of K§: 

Lemma 1 Let <3? := (ipi, (fd) G £ c , and x := (xi, ■~,Xd) G (R d ) d , then 
x G K§ if and only if for all i G {1, ...,d}: 

Y t <P j (x j ) = <p*((-l) i + 1 j\x j ) (22) 

and 

f\ Xj €d<pi(xi). (23) 

Proof. If x G Kq, then 

which proves (|22p . Now let j/j G K d , since x G -RT<j> and $ G £ c , we have: 

<ft(m) - 9ii x %) > \Vi- x u (-i) i+1 A x i / ' 

\ jy< / 



which proves (|23p . 

Conversely assume that x satisfies (j22j) - (j23j) for some i. From (j22j) . we 

get 

J^kfaO = ^*((-l) m /\xj) + <pi(xi) 

k=l j^i 

with (|23|) . this yields 

X] <Pk(xk) = ( Xi, (-l) i+1 A x i) ) = det(x), 
fc=i \ ^» / 

so that x G q 
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Remark that in fact, for each fixed i, K§ consists of those (x G W d ) d that 
satisfy f|22|) - ()23| ) for that particular index i. We then immediately deduce 
the following: 

Proposition 3 Let (ifi, ...,%) G £ c , the following assertions are equivalent: 

1. ((pi,...,<pd) solves (£>), 

2. there exists 7 G II(/ii, ...,Hd) such that for every i G {1, ...,d} and for 
7-a.e. (xi,....,x d ) G (R d ) d ; (df and ftoW, 

3. there exists 7 G II(/ii, /i<j) and i G {l,...,d} suc/i that for 7-a.e. 
(xi,....,a;d) and hold, 

If the marginals have additional regularity, we immediately deduce the 
following uniqueness result for the dual problem: 

Proposition 4 If for every i G {l,...,d}, fa is absolutely continuous with 
respect to C d and has a positive density with respect to C d then (T>) admits 
a unique solution (up to the addition of constants summing to to each 
potential). 

Proof. Let 7 be a solution of (M.K,) and $ := (cpi, (fa) solve (£>), it 
follows from the duality relation that the support of 7 is included in K§. 
Moreover, it can be assumed that each tpi satisfies (|19p hence $ G £ c . Let 
i G {l,...,d}, and let us desintegrate 7 with respect to its i-th marginal : 
7 = fa ® j Xi . We deduce from lemma [1] that for almost every X{ G M rf and 
j Xi almost every (xj)j^i G (M a! ) d ~ 1 one has: 

f\ Xj ed<pi(xi) 

hence, by convexity 

(-i) m f ^x^dx^ed^). 

Since <pi is differentiable fa-a.e., we get that for fa-a.e. Xi, one has 
Vipi(xi)= ^Xjd^dxj)^) 

since the rightmost member of this identity does not depend on <£, we are 
done. 1 — 1 
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To sum up, getting back from (£>) to {M.K), we obtain the following 
characterization of optimal transportation plans: 

Theorem 3 Lei 7 £ II(/xi, 7 solves (MJC) if and only if there exists 
l.s.c. convex functions ipi : R. d — > R U {+00} smc/i i/iai /or a// i £ {1, d}: 

> A^) °» ( Rd ) d ' ( 24 ) 

J>,(x,) = ((-l) i+1 /\s,-) 7-o-e., (25) 

(-l) i+1 /\ije%(i i ) 7 -o.e.. (26) 

Once again, one can replace "for all i € {l,...,cT}" by "for some i £ 
{1, ...,d}" and "7-a.e." by "on the support of 7" in the previous result. Of 
course, any (<^i,...,<^) satisfying the previous statements is a solution of 
(V). 

To illustrate the previous considerations, let us consider the case d = 3 
and assume that (ip, ip, x) is a (known) triple of convex potentials that solve 
(V). For the sake of simplicity, also assume that the marginals (^1,^2,^3) 
are absolutely continuous with respect to £ 3 . Then any optimal transport 
plan 7 is characterized by the extremality condition: 

tp(x) + V'(y) + x( z ) = det(x, y, z) on the support of 7. 

This implies that, 7-a.e., one has 

4>(y) + x(z) = ^*(yAz) 
(f(x) + x( z ) = ip*(—x A z) 
ip(x) +i/)(y) = x*(xAy) 

and 

V(p(x) = y A z 
< VTp(y) = -xAz (27) 
, Vx(z) = x Ay 
Note that in particular, one has: 

(x,V<p(x)) = (y,Vi/>(v)) = («>Vx(2)) 

= det(x,y,z) = det(V<^(x), V-0(y), Vx(z)) 7-a.e. . 

The previous conditions clearly impose important geometric restrictions on 
7. It implies in particular that for /xi almost every x, the conditional proba- 
bility of y or z given x is supported by V^x)- 1 . The conditional probability 
of z given (x, y) is even more constrained, indeed if (x, V<^(x)) 7^ then the 
previous conditions impose: 

Vip(x) A Vip(y) 
(x,V(p{x)) 
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4 The radial case 



In this section, we focus on (MIC) in the case where the measures (/Xj) are 
radially symmetric. This means that for all 1 < i < d, and for all R in the 
orthogonal group of W 1 , 

Rft/J-i = IH- 

Then, introducing the measures f£ = | ■ |(J/ij, we have, for all / S C c (M. d ), 

f(x)d^(x) = J ^ + °° f(re)dti(r)\ dH^e). (28) 

In order to solve (MIC), let us first remark that it is intuitive that the poten- 
tials which solve the dual problem (T>) are radially symmetric. Then, notic- 
ing that (|2*7|) generalizes to every dimension d, it follows that the support of 
an extremal measure will be included in the set of orthogonal systems. The 
only unknown here will be the relations between the norm of each vector. 
These relations will be obtained solving the following problem. 




sup ( / ( \ \r i }dY(n,...,r d )}. (MIC' 
7<- e n04\...,/^) 



Proposition 5 Assume that, for all 1 < i < d, ix[ has no atom. Then, the 
problem (MKJ) admits a unique solution given by 

dj r (n, ...,r d ) = dn\(n) <g> 5{ r2 =H 2 ( ri )} ® . . . ® S{r d =H d ( n )}, 

where, for each 2 < i < d, Hi is the monotone rearrangement map of \i\ 
into 

Proof. Proceeding exactly as for (MIC), we get the existence of a mini- 
mizer for (MKJ), together with the existence of a solution for its dual 
problem 

/ d 



™ £ „ **** > v 



{il> lt ...,tli d )ee r \f^Jo 

where £ r is the set of (tpi, . . . , ipd) S (Co(^+)) d ; sucht that, for all (r\, . . . , r d ) 
in (R + ) d , 

d 

Ti 

i=l i=l 

As for (MIC), if (-01, • • • , ipd) is such a solution, we can assume that the tpi's 
are l.s.c. convex functions on M + , and optimality conditions read as 

.fen^,...^). 
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• fory-a.e. (ri,... ,r d ) e supp^), 

Vl<i<d, ,^(r 4 )=n r i- ( 29 ) 

Let us introduce the functions Gi(r) = rip'^r). Then, multiplying (f29l) by 
rj, we get that, for 7 r -a.e. (ri, . . . , r^) in supp^), and for all 1 < i < d, 

G l {r l ) = Gi(n). 

Thanks to (|29p again, it appears that each ■0- is non negative and then, the 
function G, is non decreasing. It follows that, for ^-a.e. ri, and for all 

Vl<i<d, n = fli(n) := (Gr 1 o Gi)(n), 

where G^" 1 stands for the generalized inverse of Gj, and we get that the 
support of an optimal measure is necessarily in the closure of the graph of 
(Hi)2<i<d, where all the H^s are non decreasing. But, the constraint on 
supp(7 r ) means exactly that, for all 2 < i < d, Hi pushes fi\ forward to jj^. 
This ends the proof since, the measure \x\ being non atomic, such a map is 
unique. q 

The main result of the section is the following. 

Theorem 4 Let {ni)i<i<d be radially symmetric measures on M. d , and as- 
sume that the measures /Xj on M + have no atom. Define the measure 7 on 
(R d ) d by 



7(x) = ^i(xi) <g> 7* 1 (x 2 ) "y Xl,X2 (x 3 ) ® . . . ® 7fi>~>*dr-i 



with 



dH d ~ i ~ 1 

VI < i < d - 2, 7^'-^ = l^ (ki|)s(xij ^. ) __, 
S(x 1 ,...,x i ) = {x 1 ,...,x i } ± nS d ~ 1 , ( 3 °) 

with the maps {Hi)2<i<d given by Proposition Then, 7 is a solution to 
{MIC). 

Remark 3. The previous solution is completely explicit since the H^s are. 
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Remark 4. As in section [21 the previous construction admits a very simple 
probabilistic interpretation. Indeed, the measure 7 is the law of a vector 
(Xi, . . . , X d ) of random vectors, such that, for all 1 < i < d, [i% is the 
law of Xi. Then, (|30p means that the conditional law of Xi given X\ = 
xi, . . . , Xi_\ = Xi-i is uniform on H(\xi\)§(xi, . . . , Xj_i) for 1 < % < d — 1, 
and X d is given by 

d-l 

X d = H d (\X 1 \) /\ 
i=i 




Proof. It is sufficient to prove that 7 G n(/ii, . . . , fj> d ) and that there exists 
(tpi, . . . , tp d ) € £ such that ([2Tj) holds for all (xi, . . . , x^) in supp(7~). Notice 
that supp(7~) consists of the (x\, . . . , x d ) such that 

(xi, . . . , x d ) is an orthogonal basis of M d , , , 

for /ix a.e. x\, V2 < i < d, \xi\ = Hi(\x\\). 

Let us set 

VI < z < cZ, (^(x) = V>i(|x|), 

where (^1, . . . , ^d) is a solution to the radial dual problem (T> r ). Then, for 
all (xt,...,x d ) in (R d ) d , 

d d d 

^2ip(xi) = ^V(M) > Y\ M ^ det(xi, . . .,x d ). 

i=l i=l i=l 

In addition, If (xi, . . . , x d ) belongs to the support of supp(7~), then by (f3T1) . 
we both have 



det(xi, . ..,x d ) = Yl 



I 

i=i 



and (\xi\)i<i< d belongs to the support of the solution 7 1 " of (M.K7 ) by propo- 
sition [5j Since (ipi)i<i< d is optimal in (V r ), it follows that 



d 

{Xi) 

i=l i=l 



det(xi, ...,x d ) = ^2^i(\xi\) = ^2<fi(a 



To end the proof, we just need to prove that 7 has its marginals the /Vs. 
The first marginal of 7 is obviously Let / € C^M 6 '), then by definition 
of 7, 



dH d ~\x 2 ) \ 
H 3 (|*i|)S(xi) H 2 {\xi\) d 2 \E d 2 \ J 
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Using the radially symmetry of fix, performing the change of variables %i 
H2{\x\\)y in the inner integral, and then applying Lemma[2l we get 




/(g 2 (r)y) ^_y )d, l [ [ r) 

£iW <, - 1 (<7)\ 



by the definition of H2. Repeated applications of this argument lead to the 
same result for the other marginals of 7 up to the (d — l)-th. For the last 
one, we proceed as in the proof of Proposition [TJ For a function / S Co(M. d ), 
by change of variables, 

L. l{ rn^ d - 2 ) f { Hd(r) h (Si)) 

-L .vHSfe))^ 



f(H d (r)x d ) , 

(x 1 ,...,x d _ 2 ) P I 



noticing that, if xi, . . . , 2^-2 are fixed, the map 

d-l 



is a rotation with angle tt/2 on . . . , cc<2-2)- We can then repeat the 

argument used for the other marginals and then finish the proof. 



5 Volume maximization 

So far, we have restricted our attention to the case where the objective 
function is the determinant although we were initially motivated with a 
volume maximization problem which corresponds to: 

(MK,) a sup / \<M,(x 1 ,...,x d )\d'y(xi,...,x d ). 

76lI(^i,...,/Lt d ) JR d * d 
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instead of (MIC). In the radial case, treated in the previous section, we have 
seen that optimal measures for (MIC) give full mass to the set of matrices 
with nonnegative determinant. In this case, there is no loss of generality in 
replacing (MJC) a by (MIC) in the sense that solutions of (MIC) also solve 
(MIC) a - This holds true under less restrictive symmetry assumptions on the 
marginals: 

Proposition 6 Assume that (T7p is satisfied and that among the marginals 
(fii, fid), two are symmetric (i.e. (— id)(j/ij = /ij, (— id)jj/ij = \ij for 
two different indices i and j) then any solution of 7 of (MIC) satisfies 
det(x\, ...,Xd) > — j-a.e. and actually also solves (MIC) a . 

Proof. Without loss of generality, assume that /xi and \i2 are symmetric 
and that d > 3. Assume that (931, ifd) € £ c solves the dual problem (V), 
then so does (ip\ o (—id), if 2 o (—id), 933, (p,i). Hence (ip±, ■■■,ipd) also solves 
(T>) where ip\ and ip2 are the even part of ipi and ^2 respectively and ipi = cpi 
for i > 3. Since -01 is even we actually have for all (x\, ...,Xd)- 

d d 

ipi(-xi) +y^^j(xj) = y^ipi(xj) > I det(a;i, ...,x d )\ (32) 

j=2 i=l 

but, from the duality relations, one deduces that 7-a.e. 

d 

y~]ipi(xi) = det(xi, ...,x d ) > 0. 
i=i 

To show that 7 solves (M)C) a , we use a classical duality argument. Indeed, 
let 7 £ n(/ix, ■■■,(J,d), from (f32l) . one deduces: 

/ I det(xi, ...,x d )\dj(xi, ...,x d ) < ^^{(x^d^Xj) 

jRdxd 




det(xi, ...,x d )\drf(x!, ...,x d ). 



□ 



6 Uniqueness issues 

We end the paper by some remarks on uniqueness. For this, we will mainly 
use the very simple example of section [2j First, as noticed at the end of 
section [21 if the objective function is 

x 1— ► I det(x)|, 
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then there are infinitely many possible choices for the third random vector, 
given the first 2 ones. Indeed, for all x and y in the ball B with \x\ = |y|, 
and for all p £ [0, 1], then we can choose for the conditional law of Z, given 
X = x and Y = y: 

p5xAy + (1 —p)5yAx. 

M R 

Note that for the determinant, only the value p = 1 is optimal and it cor- 
responds to the solution given by Theorem [TJ In the sequel, we shall work 
with the determinant as objective function. 

As seen in the introduction, in dimension 2, since the problem (MIC) can 
be reduced to the classical Monge-Kantorovich problem with quadratic cost, 
the solution is unique, at least for sufficiently regular measures. Actually, 
we have more, that is, up to a rotation, the optimal measure is supported by 
the graph of the gradient of a convex function cp. In terms of random pairs 
(X, Y), the optimal coupling is characterized by the relation RY = Vip(X), 
with (p convex and R the rotation with angle —tt/2. In higher dimensions, 
such uniqueness is lost. 

First, let us consider the solution given by theorem [TJ For x E B, and 
given X = x, we define the conditional law of Y by the measure on |x|S(a;), 
absolutey continuous with respect to the 1-dimensional Hausdorf measure, 
with density 

I V I 

where e is an arbitrary fixed vector in S 2 , which is a probability measure 
since we add to 1 a function which is greater than —1 and has zero mean on 
any circle S(x). The vector Z is then chosen exactly as in Theorem [TJ The 
proof of the optimality of this measure is the same as for Theorem [TJ 

Finally, we study the possibility for the extremal measure to have the 
same structure as in the 2-dimensional case, that is to be supported by 
the graph of a map T : M. d — > Such a T would then be called a 

"Monge" solution. Notice that the solution given by Theorems [U and [T] does 
not have this structure but the question of the existence of such a solution is 
quite natural. Let us set the problem on M. d with all the measures uniform 
on the ball, which is actually covered by Theorem [TJ Then, according to 
the extremality conditions, the existence of a Monge solution for (MIC) is 
equivalent to the existence of (d — 1) maps T2, ■ ■ ■ T^-i- B — > B such that: 

for a.e. i£8, {x, T%(x), . . . , Td(x)} is an orthogonal basis of M. d , 
for a.e. x £ B,\/2 < i < d, \Ti(x)\ = \x\, 
V2 < i < d, Ti preserves the Lebesgue measure. 

(33) 

The construction of such maps seems to be very difficult in general and 
we left it open. However, in dimension 4, we have the following explicit 
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construction: 





-X2 


\ 




( 


-X3 \ 




X\ 




, T 3 (x) = 




X4 




— X4 






X\ 


\ 


X3 


J 




K 


-X2 J 



T 4 (x) 



\ 



-X4 \ 

-X3 

X2 

X\ ) 



In contrast, in dimension 3, it is not clear to us whether one can even find 
a single measure preserving map T : B — > B such T is norm-preserving and 
a.e. orthogonal to the identity map. If such a map exists, then by the hairy 
ball Theorem, it cannot be continuous. 



Appendix A 

Lemma 2 Let k > 2 be an integer and f G Co(S k x S k ). Then, 

dH k -Hy) \ dH k {x) _ f ( f dH k -Hx) \ dH k (y) 

(34) 



ls k US(i) I s I / I s I J§ fc \JS(y) 

where 



Proof. Let us introduce the following notations. Let us define the measure 
fi on S fc x § k by 

dfi(x, y) = dv(x) ® d6 x (y), 

where 

, . . (iW fc (x) , ,dH k -\y) 

dv{x) = v and = l g(3 .)(y) igfc-ii • 

The identity (j34"|) expresses the fact that is symmetric in x and y. We set, 
for 2; € S fc and e > 0, 

§ e (x) = |z G S fc ; I (z,x) I <e} . 

Let us notice first that 7^ fc (S £ (x)) is independant of x and denote this quan- 
tity by u e . Indeed, for all x and x' in S fc , let R be a rotation that maps x 
to x'. Then, we have 

S £ (x') = S £ (i?x) = R- 1 (S £ (x)) , 

and then Tt k (§i £ (x')) = H fc (8 £ (x)) since the Hausdorff measure is invariant 
under rotations. Next, we define the probability measure 6 X on § fc as 

^(y) = %(. ) (y)^ M - 

We 
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For all / e C (S fe ,S fc ), and for all x e S k 



[ f(x,y)df)*(y)- I f(x,y)dO x (y) 
Js k Js k 



< 



(f(x,y)-f(x,T x (y)))de x (y) 



+ 



f(x,T x (y))d8 x (y)- / f(x,y)d0 x (y) 
Js k 



(35) 



where 



T x {y) 



Wxy\ 



and 7r x is the orthogonal projection onto x^. First, we shall estimate the 
first term in the right-hand side of (j35j) . It is not difficult to see that if 
e < 1/2, then for all x <E S k and y G S £ (x), 

\\y-T*(y)\\<4e, 

thus, denoting by ujf the modulus of continuity of /, 



(f(x,y)-f(x,T x (y)))d0 x E (y) 



< LUf(4s) 



For the second term in the right-hand side of (j35[) . we write 



f{x,T x {y))d9 x {y) = / f(x,z)d(T*m(z), 
s k Js k 



where 



T x (y) 



Since T x commutes with any rotation that leaves x invariant and 9 X is 
invariant under such a rotation, the measure T X %Q X has its support in §(x) 
and is invariant under rotation on §(x). So necessarily, we have 



and 



T xpx = gx 



f(x,y)d0 x (y)- / f(x,y)d6 x (y) 
s k Js k 



< u f (4e). 



As a consequence, 



lim sup 

\x& k 



f(x,y)d9 x (y)- / f(x,y)d6 x (y) 
S k Js k 



0. 



(36) 



We set dfi £ (x,y) = dv(x) ®d9 x {y). Then, (f3"6"|) implies that (n e ) e converges 
weakly to /i. In addition, fi £ is symmetric in x and y, which ends the proof. 



□ 
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Appendix B 



Proof of theorem [2j 

Stepl : existence of maximizers for (MK) 

Let fceN and define: 

H k := max(H, -k) (37) 

By construction, —k< H_ k < and H_ k decreases to H . 

Let 7 n be a maximizing sequence of (MIC). Since H(ni, ■ fJ>d) is tight, it 
follows from Prohorov's Theorem, that (taking a subsequence if necessary) 
there exists 7 € M(R dxd ) such that for every continuous bounded function 
/ on R dxd , J Rdyd /7 n converges to j Rdxd fdy. Clearly, 7 G II(/zi, and 
for every k 

H k d^ = lim / iLk^ln — hm / Hdj n = sup(MJC). 

id n J R dxd n J R dxd 

It thus follows from the monotone convergence theorem that 7 is a solution 
of (MIC). 

Step 2 : duality by approximation 

Consider the problem: 

(MlC k ) sup / H k (x)dj(x) 
yeTl(pi,...,iJ, d ) JR dxd 

By the same considerations as in step 1, (MJCk) admit solutions, let 7^ be 
such a solution. Now for e >, let r = r £ > be such that 

d „ 

S2fi t (R d \B d )+ / (l + \H (x)\)d-f(x) < £ , V 7 6n( W ,.., W ). (38) 



Define then for i = 1, d: 



Consider now: 



(MK. T)k ) sup / H k (x)d , y(x) 



7gn(^t 1 ' ,...,/i d - ) ■ 

and (defining £ r ,fc as in (fT2|) with H_ k instead of det) its dual: 



d „ 

(2? r>fc ) inf V / ^d/i[' fc . 
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We know from proposition [2] that the values of (MtC r ^) and (XV.fc) are 
equal and both attained. Moreover, using again proposition [21 the infimum 
in (XV,fc) is attained by some ip r,k = ((p r ^ k , (f d k ) such that: 

<Pi' k (vi)= SU P < £Lk(xi, ...,x d ) -^^(xj) > , Vxj £ B r , (39) 

and, for all x £ B r : 

-k< <fi k (x) < 0, < $ k {x) <k, i = 2, d. (40) 
Define then for all X\ € R rf : 

ijj r { k ( Xl ) = sup I E h (xi, x d ) - ^2 Vj k (xj) } ( 41 ) 

(x 2 ,....,x d )aBi- 1 y j=2 J 

Construct then inductively ify ; ■••> V'S-i by setting for i = 2, — 1 and 

xi e R d , 



^ k {xi) = swp{H k (xx, ...,x d ) - ^2^ r j' k (xj) - ^tf 

j<i j>i 

(x 1 ,...,x i . 1 ,x i+1 ,...,x d ) eR dx ^ xB d ^} 



(X 



Jh 



and finally 



V^Od) = SU P 



H k (xi,...,x d ) - ^^ k {xj) 

3=1 



It can be checked easily that each ip^' k extends to the whole of M. d . 
By (j4T)1) it can also be checked that ||V' r ' fe ||oo = 0{k). By construction 
ifj r ' k = (i/}^ ,k , ... : (p d ,k ) E £ fc , where £ k is the set of d-uples of lower-semi 
continuous functions {ip%, <p d ) from R d to R U {+00} such that 

d 

^2v>i(xi)>Hj,(xi f ...,x d ), y( Xl ,...,x d )eR dxd . (42) 

i=i 

In other words, £ k is the admissible set of the dual problem of (A41C k ): 

d „ 

(V k ) inf V" / ipidm. 
Let j r>k be a solution of (AiK r> k) and define: 

7r,fc : = lk(B d )~f r ,k + l K dxd\ B d. 
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It is easy to verify that 7 rfc 6 ■■■,fJ*d) hence we get: 

sup(.M/C fc ) = / H_k d lk > / Kkdlr,k 

=lk{Bf.) [ H k d lr}k + [ H k d lk 

JBi JtL dxd \Bi 



-~lk{ 



\B d r 

Using (j38l) . we thus get: 

sMMKk) > lk {B d r ) ( E / V? V* ) " e = E / 

\i=i / i=i 

d „ d „ 

= E / v^Tfc - E / - e > inf (^) + 

Letting e go to we then get sup(A4/Cfc) > inf(Pfc). Using Prohorov's 
Theorem, taking a subsequence if necessary, we may assume that there exists 
7 6 .., /Xrf) such that for every continuous bounded function / on M dxd , 

j R dxd flk converges to j R dxd fdy. Since H_ k > H_ = det — Hq, we then have 
for all k: 

sup(MICk) = / Kkd-Jk > inf(P) - a 

J^dxd 

Since H_ k is nonincreasing in k, for all ko, we then have: 

/ E_k dl = lip / Kk dlk > / Kk d lk > inf(D) - a 

J^dxd k J^dxd JRdXd 

taking the infimum in ko and using the monotone convergence theorem we 
then get: 

sup(MlC) -a> f Hdry > inf(P) - a. 

Since we already know that sup(7W/C) < inf(P), we have sup(.M/C) = 
inf(Z>). 

Step 3: existence of minimizers for (P) 

It remains to prove that the infimum is attained in (T>). By the convexi- 
fication trick, we can find a minimizing sequence of (T> r ), ip n := (tp™, 
such that for all n, all i and all X{ G one has: 



d \ d 



<Pi(xi) = sup 

(x j ) j9 , i m dxd - 1 



det (zi,..., s,0 - E ?i foOf ( 43 ) 
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Note the functions iff are convex l.s.c. and not identically equal to +00. 
By (03]), (pf admits therefore some affine minorant hence ^{x) + \x\ Pi /pi 
achieves its minimum. Using again that the objective in (£>) is unchanged 
when changing (ipi, ip d ) into (ip± + a%, (fd + a d ) for constants ai that 
sum to 0, we may then also assume that for all i = 1, d — 1: 

min lifi{x) + — 1 = 0. (44) 



We easily deduce from the previous and (j4*3|) that 

\x\ Vd 

<Pd(x) + U — > 0, Vx G R d . 
Pd 

Let m > and define truncated potentials </?™' m and V>"' m by: 



C"» = ^ ,m (x) + ^- = min((^™(x) + m) Vx G R d . (45) 



,n,m/ \ n,m/ 

' x ; -r - 

Pi Pi 



By construction, < ^ ' < m and there exists C such that for all (n,m), 
one has: 

£ / C' m ^ < C. (46) 
and since V>"' m > 0, for all (xi, ...,Xd) G IR a!xc( , one has: 

d 

y^' m {xj) > min(#(xi, ...,x d ),m). (47) 
i=i 

For fixed to, taking subsequence if necessary, we may assume that tpi' m 
weakly converges in L 1 (/ij) to some limit ip™ as 11 — > +00. By Mazur's 
Lemma, there is some sequence of convex combinations of the (ip l i ' m )i> n that 
converges in L 1 ^) and (possibly after an extraction) /^-almost everywhere 
to as n — > +00, there is no loss of generality of setting ip™ = +00 outside 
the set where there is convergence. Hence, passing to the limit in (|47p . we 
get, for all (xi,...,x d ) G R dxd 

d 

Y,^T{xi) > rmn(H( Xl ,...,x d ),m). (48) 
i=i 

Now, ip™ is a nondecreasing sequence of nonnegative functions that is bounded 
in L 1 (//j). By the monotone convergence theorem, ipf 1 converges pointwise 
and in L 1 ^) to some function ^j. Defining ifi{x) := ipi(x) — \x\ Pl /pi, and 
passing to the limit in (jlBl) first yields: 



5>™(^) >det( Xl ,...,x d ), V( Xl ,...,x d ) €R dxd . (49) 



i=i 
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Since ipf > ^ m ', we also get: 



d „ / d ,. 

inf(P) =lim V / ^d/^ > sup lim V / 

d „ d „ 

= ^idm - a = y~] / <^id//i 



Thanks to (|49l) . we can apply the convexification trick to ((pi, ...,<pd). We 
then obtain an element of £ that satisfies (|43p and solves (P). 



□ 
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